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We consider the stability of three Coulomb charges {+1,-1,-1} with finite masses in the 
framework of nonrelativistic quantum mechanics. A simple physical condition on masses 
is derived to guarantee the absence of bound states below the dissociation thresholds. In 
particular this proves that certain negative muonic ions are unstable, thus extending the old 
result of Thirringi to the actual values of all masses. The proof is done by reducing the 



initial problem to the question of binding of one particle in some effective potential. 
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INTRODUCTION 

> ' 

The stability of three particles with pure Coulomb forces is an old and extensively studied 
problem, this is to explain why certain ions and molecules stay as a whole and some dissociate into 
a bound pair and a single particle. Under stability of the Hamiltonian H we shall understand the 
existence of a bound state with the energy strictly less than inf a ess (H ) , i.e. a stationary state 
below all dissociation thresholds. As a quantum system three particles with Coulomb interactions 
demonstrate interesting behavior. It is known that three charges {1 + e, — 1, — 1} for any e > 
form the system, which is stable regardless of mass values. However at e = the situation abruptly 
changes and due to the screening effect not all systems remain stable. The typical example is the 
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unstable muonic hydrogen ion p/s e , where the heavy muon is tightly bound to the proton and 
screens the positive charge for the electron. This interesting effect is well studied, in particular it 
has been proved 2 that for equal dissociation thresholds the system remains 
very well how stability depends on masses and charges). 

Among the instability proofs for three unit charges the most appealing is that of Thirring 1 , 
which does not require any numerical calculations. Thirring considered a negative hydrogen ion 
with an infinitely heavy proton, and proved that such negative ion is unstable when its second 
negatively charged particle is lighter than electron by a factor larger than it (tt can be replaced^ 
with a better constant 1.57). In his method Thirring exploited the fact that the ground state in 
the hydrogen atom is substantially separated from other states in its spectrum (the non-degenerate 
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energy levels have 1/n 2 dependence), and thus its role becomes emphasized. This suggests an idea 
to move the problem to the space spanned by the projector Pq = \(f>o){(f>o\ 8> 1, where 4>o is the 
ground state of Hydrogen. After estimating the part of the repulsion, which is present in this 
space, the problem reduces to checking the binding of one particle in some effective potential. The 
contribution from the attractive interaction term coming from additional particle is easy to treat 
because it commutes with Pq. Yet this is no longer true if one considers particles of finite mass or 
any system of four particles. 

Thirring's bound was improved 3 and it was shown that the muonic ion pfi~~e~" is unstable for 
actual values of all masses. However this extension (Eq. 24 in Ref. 0) is weak in the sense that 
it fails when the second particle is heavy compared to other particles. (For a physical example, 
it does not prove that the ion [i~pe + is unstable, yet we shall prove it here). This extension still 
uses Thirring's treatment of repulsion and it is unclear how one could extend it to four particles. 
Armour 5 with his method proved the instability of such systems as positron-hydrogen-atom e~pe + 
and fi~pe + , but this method requires certain numerical assistance. It also bases on the separation of 
variables in the problem of two fixed centers, which makes it inapplicable to four particles because 
the variables do not separate even in the case of three fixed centers. 

Here we follow the Thirring's idea but the nucleus does not have to be infinitely heavy. The 
derived physical condition restricts the ratio of Jacobi masses, which makes the system stable. It 
can be used in conjunction with Thirring's result and convexity properties of stability curve 3 for 
a reasonable determination of stability area free from any numerics. The present method has an 
important advantage in that it admits generalization to four particles. It should be mentioned that 
both Thirring's method and the new one share the same deficiency, namely, they are not applicable 
when the dissociation thresholds in the system are close or equal (e.g. when two like charges have 
equal masses). In particular, using such methods one cannot prove the "overheating" effect, when 
a system with charges {+1, — q, —q} and equal masses loses its binding for some q > 1. Nowadays 
various charged particles are produced in laboratories and it is, of course, of interest to know the 
principles behind formation of exotic atoms or molecules. All this motivates the present analysis. 

II. FROM THREE PARTICLES TO ONE IN EFFECTIVE POTENTIAL 

Let rrii,qi,Ti 6 M 3 denote masses, charges and position vectors of particles i = 1,2,3. We 
put q\ = +1, and </2,3 = — 1, and the interactions between the particles are Vik = %%/|rj — r&|. 
We enumerate the particles in such a way, that the particles (1,2) form the lowest dissociation 
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threshold. The stability problem with Coulomb interactions is invariant with respect to scaling 
all masses^, so we can put H = 1. We separate the center of mass motion in the Jacobi framed 
putting x = ri — r2 and y = r3 — ri + ax, where a = m%l{m\ + 7712) is the mass parameter 
invariant with respect to mass scaling. The reduced masses and Jacobi momenta are respectively 
fJ-x = mimi/ (jn\ + 1712), li y = 7713(772-1 + 1712) /(mi + 777,2 + m 3 ) and p x>y = —iV x ,y It is convenient 
to scale all masses so that fj, x = 2. (In Sec. [El we shall rescale them back). The Hamiltonian for 
the system on the tensor product space L2(R 3 ) <8> L2(M 3 ) is 

H = hi2®l + l®P L + W 

2\Xy 

where 



W = Vl 3 + V 23 = — r + r — -r (1) 

|ax-y| |(l-a)x + y| 

and hi2 = P x /^ — 1/x is the Hamiltonian of the pair of particles (1,2) (notation x is used instead 
of |x|). The ground state wave function of h\2 is 00 = \/8/tt exp(— 2x) so that /ii20o = ~ ( t ) o- 
The Hamiltonian H is self-adjoint on T>(H) = H 2 (M. & ) (square integrable functions having partial 
derivatives up to the second order in the weak distributional sense) and by the HVZ theorem 
a ess (H) = [—1, 00). We split positive and negative parts of W by introducing W- := (\W\ — W)/2 
and W + := (-W)- = (\W\ + W)/2 and we have W = W+ - W-, where W± > 0. Instead of H we 
shall consider the Hamiltonian 

77 2 

H = hi 2 ® 1 + 1 <g> -p- - W. (2) 

Note that the part W- can also be expressed as 

W- = -(Vi 3 + V 23 F(x,y)) (3) 

where F = 1 when | V13 1 > | V23 1 and F = | V^L3 1 / 1 V23 1 when | V13 1 < | V23 1 , and we have < F < 1. 
Because ||V23-F</>|| < ||^30|| we can directly apply Kato's theorem^ on self-adjointness of atomic 
Hamiltonians and find out that H is self-adjoint on T>(H) = D(H). We cannot though directly 
apply the HVZ theorem to locate inf a ess (H) but we observe that the following inequality holds 
(H — V23) < H < H . Here (H — V23) is the original Hamiltonian without repulsion, which is 
bounded from below and by the HVZ theorem va.ia ess (H — V23) = — 1. Thus from the min-max 
principle^ we get inf a ess (H) — — 1. 

Now let us assume that H is stable, i.e. H has a bound state with the energy less than — 1. 
Because H < H from the variational principle we conclude that H also has a bound state \P G D(H) 
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with the energy below inf a ess (H) = — 1 which means 

= (-l-8)V (4) 

where 5 > is the extra binding energy. Let us introduce a projection operator Pq = |</>o)(</>o| ® 1 
(P : V(H) -> £>(#)) and put n := P ^> and £ := (1 - P )^, where obviously r/_L£ and ^ = r/ + £ 
and r/, £ E T>(H). Taking the scalar product of each side of @ with r? and £ we obtain 

2 

® If - ^-1^ " fol w -IO = -*IMI 2 ( 5 ) 

Zfiy 

2 

m 12 ® iio + <eii ® ^ - - <eiw_i^> = (-1 - 5)na 2 (6) 

where we have used (rj\l (g) = because Pq and 1 <8> Py commute. We shall assume that 

IMIi ll£ll ( we shall get rid of this assumption in Theorem then we are free to choose such 
normalization of that ||£|| = 1. 

From the bound spectrum of h\2 we have— h\2 ® 1 > — -Po ~~ 1/4(1 — Pq), hence for the first 
term in © we get the bound (£|/ii2 ® 1|0 — —1/4. Introducing two non-negative constants 
a := \J (r]\W- \rj) and /? := ^/ (£|W_|£) we get by virtue of the Schwarz inequality |(^|W_|^)| < a/3. 
Now we can rewrite Eq. ©~® 1° obtain the main pair of inequalities 

<r?|l <g> A. - W_|fj) - a/3 < (7) 

2 

<fli®^-w_|0-a/j < -- (8) 

Using the second inequality we shall find max/3/a and by substituting this value into Q we shall 
formulate the stability condition. 

Lemma 1. Suppose that Eq. Jj^ ZioWs and //„ < 3/2, i/ien i/ie following inequality is true 

P<( x [JL-i) a (9) 



2 fx 
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Proof. First, let us show that for A > 

inf Ml AW-\ X ) > (10) 

llxll=i 

It suffices to prove this for x £ C^°(R 6 ). Using Q from the variational principle we get 

A \ , x 



2 / 2 

(x|i®^-aw_|x> ^/ dx / ^yx*(x,y) I ^ 



> 



|ax-y| 

A2 ^| dx J dy| X | 2 (x,ax + y ) = -^|| X || 2 (11) 



Xlx, y) 
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from which ()1U|) follows and where we have used the explicit expression for the ground state energy 
of the hydrogen atom. (Using an appropriate set of trial functions it is easy to show that there is 
an equality sign in but we do not need this for our purposes). Now using (|10|) we obtain the 
following chain of inequalities 



2 2 

inf (£\l®p--W-\0 = max inf <g> p- - (A + + A/? 2 ) 

t£D(H) 2fi y ' a>-i £ev(H) 2fJLy 

ll?ll=i ' 11411=1 

(£,\W-\t)=P 2 <£|iy_|0=/3 2 



11=1 



/3 4 



/? 2 



2^ 

Substituting this result into © and putting a = s0 we obtain 

Now simply minimizing the left-hand side of Q12JI over 2 we obtain the lower bound on s, which 
gives us Eq. ©. □ 

It is worth noting that the relation in Lemma is a version of the uncertainty principle, when 
2 grows, the kinetic energy term grows faster like /3 4 . Let us introduce an effective potential 
V e ff(y) := f dx\(j)o\ 2 W-. We formulate the result as 

Theorem 1. If the system of three charges is stable and fly < 3/2 then the particle with mass ji y 
must have a bound state in the potential — (l + (y/3/2fi y — V e ff. 

Proof. We have 7^ 0- The function rj has the factorized form rj = 0o(x)/(y), where 

/ G H 2 (M. 3 ), Il/H / 0. By substituting © into © and using expressions for a 2 and r] we get the 
necessary condition for stability 

^y ' < (l3) 

In the next section we shall study V e ff(y) and show that it is a continuous function decaying like 
1/y 2 - Inequality (jl3j) means that a particle having mass \i y has a bound state in this potential. 

Now let us complete the proof considering the case when either ^ = or ?j = 0. If £ = we 
have = and substituting this into © we get a condition more stringent than Eq. (jl3fl . If 7] = 
we have ||£|| = 1 and a = 0, substituting this into © and using (flO|) for fj, y < 3/2 results in the 
contradiction. □ 
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III. BINDING IN EFFECTIVE POTENTIAL 

In this section we shall analyze the effective potential and find out at which values of the 
coupling constant A the Hamiltonian p 2 — XV e ff may have bound states. It turns out that the 
effective potential in our case has a nonphysical term, which is a long-range attraction of the type 
1/y 2 . This nonphysical behavior stems from cutting off the positive part of the potential and results 
in the infinite number of bound states at the point of binding (that is why it is meant nonphysical). 
But as it is well-known since the result of Hilbert and Courant 9 , even this long-range attraction 
does not guarantee binding, for A max,, y 2 V^ff(y) < 1/4 the inequality p 2 — W e f / > holds, i. e. no 
binding occurs, (see also the proof in Ref.H) vol- 2). Thus the non-trivial critical coupling constant 
exists, and we have to determine it. On the other hand, in such potentials short-range repulsive 
terms do not play any role for binding 10 . It would be of interest to get rid of this nonphysical 
behavior in the future (pay attention that the Thirring's effective potential 1 behaves at infinity like 
1/y 3 !). 

To calculate V e ff we must cut off the positive part of W. From (JJ) W < is equivalent to 
cos 6 > x/(ujy), where u> = (a — 1/2) _1 and cos# = x • y/xy. We shall consider separately two 
cases o > 1/2 and a < 1/2. 

A. Case when a > 1/2 

The integration is simpler in this case. After direct integration in spherical coordinates over the 
area where cos# > x/{uy) we obtain 



Jo V a (l ~ a ) a I — a I 

[ ds se- 4sy (^/a{l - a)s 2 + 1 - 1) + U (14) 
J o 



16y 2 
a(l — a) 

where 



U = 32y 2 dsse- isy {l/a- s) < (15) 

Jl/a 

Now we do not have to carry out integration in (|15|) . it is enough to see that it is a short-range 
repulsion, which does not play a role in our case. To calculate the first term in (|14|) we use 



sfa(l - a)s 2 + 1 < 1 + a(l - a)s 2 /2 to get 



V eff <8y 2 l dss*e-^<^ (16) 
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where after the integration we have dropped the short-range negative terms. Finally, we have 
Py — XV e ff >pl~ A(3/16)y -2 . The following inequality^ holds p 2 — (l/4)y~ 2 > 0. Thus in the case 
of binding, i.e. when such / exists that {f\p 2 — XV e ff\f) < 0, we must have A > 4/3. Comparing 
this with (|13f) we obtain that three charges form unstable system if fi y < 3/2 and 

2py (l + (^3/2^ - I)" 1 ) < 4/3 (17) 

Solving this simple inequality tells us that the system is unstable when /j, y < 2(11 — 2\/l0)/27 ~ 
0.3463. 



B. Case when a < 1/2 



First let us take a < 1/2. We shall write W(a) instead of W to point out the dependence 
on parameter a. We can alleviate the integration noting that W(a) = —W(l — a), thus we 
have W-{a) = (-W(l - a))_ = W + (l - a) = W(l - a) + W_(l - a). From this we conclude 
V e ff(a) = —W(a) + V e ff{\ — a). The additional integral W(a) := J dx.\(fio\ 2 W is easy to calculate 
and V e ff(l — a) for a < 1/2 we have already calculated. We obtain 



-W{a) = 16y 

-l/a 



co / / 1 

2 ' dse" 4 ^ " 



s - 



a 1 — a 



+ 



1 




1 


) 


S 1-a 




s 

a 





32y^ 



[ a ds e" 4 ^ f s 2 - -^-] + 32y 2 f°° ds e^s (-- —) 
Ji/(i-a) V I- a J J 1/a \a I- a J 



(18) 



'V(l-o) \ i- — uj Ji/ 

Using (|14j) . (fT^|) and approximation for the square root \J a(l — a)s 2 + 1 < 1 + a(l — a)s 2 /2 gives 
us 



K//(l ~a)<8y 2 [ ^ ds s 3 e~ 4sy + 32y 2 / " ds e~ 4sy ( — s 2 ] 

Jo JlHl-a) \l-a J 



(19) 



Summing ()18|) and (|19j) and calculating the integrals explicitly gives us the following expression 
for a < 1/2 



V eff < 



-Ay/a 



16y 2 



8(1 -a)y h2 2-a | 1 



-2yw(u; + 2) 2 + (3w/2 + l)(o> + 2) - -(3w/4 + 1) + 



y ' 1% 2 



(20) 



For a < 1/2 we have uo < —2, and it is easily seen that all terms in square brackets are positive 
(this leads again to short-range potentials), meaning that V e ff< (3/16)y -2 , which gives the same 
condition for stability as (|17j). We do not consider explicitly a = 1/2, it is done analogously and 
also results in (fT7|) . 
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IV. SUMMARY 

We have initially scaled all masses — ► Imij \i x , making fi x = 2. Now rescaling it back we 
get through (|17[) that the system of three charges is unstable if Py/n x < (H ~~ 2\/l0)/27 — 0.1732. 
In the case of infinitely heavy nucleus this is 7713/777-2 < 0.1732, which is worse than the refined^ 
Thirring's estimate 7773/7772 < 1/1.57. The accuracy is lost at the point of cutting the positive 
part of the potential, which induces a long-range attraction. However this is more than enough to 
prove that the muonic ions pfi~e~ or fi~pe + are unstable for the actual values of all three masses. 
The case of four unit charges {+1, +1,-1,-1} is treated similarly but the calculations are more 
involved^ and results would be published elsewhere. Let us also stress that the obtained condition 
is physical. Both Jacobi masses determine Bohr radii for the particle orbits, the orbit within the 
pair of particles (1,2) and the orbit for the third particle in the field of this pair with respect to the 
pair's center of mass. If the orbit of one negative particle is outdistanced then the attraction from 
the positive charge is screened off by the other negative particle and the system becomes unbound. 
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